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With a sample of (225.2±2.8) x lO" J/ip events registered in the BESIII detector, J/ip -5> 77r+7r"r7' 
is studied using two r;' decay modes: r;' — )■ tv'^tt't] and rj' — >■ ■yp". The X(1835), which was previously 
observed by BESIl, is confirmed with a statistical significance that is larger than 20cr. In addition, 
in the -k^-k't]' invariant mass spectrum, the X(2120) and the X(2370), are observed with statistical 
significances larger than 7.2g and 6.4(7 , respectively. For the X(1835), the angular distribution 
of the radiative photon is consistent with expectations for a pseudoscalar. 

PACS numbers: 12.39.Mk, 12.40.Yx, 13.20.Gd, 13.75. Cs 



A 7r+7r~?7' resonance, the X(1835), was observed in 
J/^ — > 77r+7r~77' decays with a statistical significance 
of 7. To- by the BESII experiment [l|. A fit to a Breit- 
Wigner function yielded a mass M = 1833.7±6.1(stat) ± 
2.7(syst) MeV/c^, a width T = 67.7 ± 20.3(stat) ± 
7.7(syst) MeV/c^, and a product branching fraction 
B{J/il^ -fX) ■ B[X TT+TT-rji) = (2.2 ± 0.4(stat) ± 
0.4(syst)) X 10^"*. The study was stimulated by the 
anomalous pp invariance mass threshold enhancement, 
that was reported in J/-0 jpp decays by the BESII 
experiment [2] and was recently confirmed in an analysis 
of V'' TT+TT" J/-0, J/tjj jpp decays by the BESIII ex- 
periment 3. The possible interpretations of the Ar(1835) 
include a pp bound state \^^, a glueball a ra- 

dial excitation of the 77' meson [llj, etc. A high statistics 
data sample collected with BESIII provides an opportu- 
nity to confirm the existence of the Ar(1835) and look for 
possible related states that decay to Tr+Tr^Ty'. 

Lattice QCD predicts that the lowest lying pseudo- 
scalar glueball meson has a mass that is around 
2.3 GeV/c^ fl?]. This pseudo-scalar glueball may have 
properties in common with the rjc , due to its similar 
decay dynamics that favor decays into gluons. One of 
the strongest decay channels of the rjc is Tr+Tr^ry' . 
Thus J/ip — >■ j'K~^TT~r]' decays may be a good channel 
for finding ^ glueballs. 

In this letter, we report a study of J/^ — >■ "fTT'^TT~ri' 
that uses two 77' decay modes, 77' — ^ and 77' tt^tt^t]. 



The analysis uses a sample of (225.2 ± 2.8) x 10^ J/^p 
events [^^ accumulated in the new Beijing Spectrome- 
ter (BESIII) jT3| located at the Beijing Electron- Positron 
CoUidcr (BEPCII) [15] at the Beijing Institute of High 
Energy Physics. 

The design peak luminosity of the BEPCII double-ring 
e+e^ collider, is 10'^'^ cm^^s^^ with beam currents of 
0.93 A. The BESIII detector has a geometrical accep- 
tance of 93% of 47r and consists of four main compo- 
nents. 1) A small-celled, helium-based main draft cham- 
ber (MDC) with 43 layers. The average single wire res- 
olution is 135 /im, and the momentum resolution for 
1 GeV/c charged particles in a 1 T magnetic field is 
0.5%. 2) An electromagnetic calorimeter (EMC) com- 
prised of 6240 Csl (Tl) crystals arranged in a cylindrical 
shape (barrel) plus two endcaps. The energy resolution 
for 1.0 GeV photons is 2.5% in the barrel and 5% in the 
endcaps, and the position resolution is 6 mm in the barrel 
and 9 mm in the endcaps. 3) A Time-Of-Flight system 
(TOE) for particle identification composed of a barrel 
part with two layers with 88 pieces of 5 cm thick, 2.4 m 
long plastic scintillators in each layer, and two endcaps 
each with 96 fan-shaped, 5 cm thick, plastic scintillators. 
The time resolution is 80 ps in the barrel, and 110 ps in 
the endcaps, corresponding to a 2a K/tt separation for 
momenta up to 1.0 GeV/c; 4) A muon chamber system 
(MUG) made of 1000 m^ of Resistive Plate Chambers 
(RPC) arranged in 9 layers in the barrel and 8 layers in 
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the endcaps and incorporated in the return iron of the su- 
perconducting magnet. The position resolution is about 
2 cm. 

Charged-particle tracks in the polar angle range 
cos 6*1 < 0.93 are reconstructed from hits in the MDC. 
Tracks that extrapolate to be within 20 cm of the inter- 
action point in the beam direction and 2 cm in the plane 
perpendicular to the beam are selected. The TOF and 
dE/dx information are combined to form particle iden- 
tification confidence levels for the n, K, and p hypothe- 
ses; each track is assigned to the particle type that cor- 
responds to the hypothesis with the highest confidence 
level. Photon candidates are required to have at least 
100 MeV of energy in the EMC regions | cos 6*1 < 0.8 and 
0.86 < |cos0| < 0.92 and be isolated from all charged 
tracks by more than 5°. In this analysis, candidate events 
are required to have four charged tracks (zero net charge) 
with at least three of the charged tracks identified as pi- 
ous. At least two photons (three photons) are required 
for the 77' — s> 7/3 (77' — > tt+tt^t]) channel. 

For J ftp — >■ ^TT'^ rj' [rj' ip), a four-constraint 
(4C) energy-momentum conservation kinematic fit is per- 



formed to the 777r 



hypothesis. For events with 



more than two photon candidates, the combination with 



is used, and 



xlc 



< 40 is required. 



the minimum x 
Events with \M^j 

0.03 GeV/c^, 0.72 GeV/c^ < M^^ < 0.82 GeY/. 
press the background from 7r'^7r"''7r~7r"'"7r~ , t/tt^ 



m^o| < 0.04 GeV/c^, \M^^ - m^l < 

-y-y \ KJ.U^ VJC V ^ U Or 

< 0.007 GeV/c^ are rejected to sup- 
and J Tr~^ '!T~ 7] [rj — >■ 77r+7r^), re- 



combinations shown in Fig. [ija). Candidate p and rj' 
mesons are reconstructed from the tt^tt^ and 



< 0.2 GeV/c2 and \M 



pairs with |M^+^- — 
rrirj'l < 0.015 GeV/c^, respectively. If more than one 
combination passes these criteria, the combination with 
M 



closest to rriri' is selected. After the above se- 
lection, the X(1835) resonance is clearly visible in the 
7r~^Tr~ri' invariant mass spectrum of Fig. [TJb). Also, ad- 
ditional peaks are evident around 2.1 and 2.4 GeV/c^ as 
well as a distinct signal for the rjc- 

For J/tp — > jTT^TT^rj' {rj' — ^ tt+tt"/;), a 4C kinematic fit 
to the 7777r+7r^7r+7r^ hypothesis is performed. If there 
are more than three photon candidates, the combination 
with the minimum xic selected, and xic < 40 is re- 
quired. In order to reduce the combinatorial background 
events from 7r° 77, \Mj^ ~ m^o\ > 0.04 GeV/c^ 
is required for all photon pairs. The rj candidates are 
selected by requiring |M^^ — m,,| < 0.03 GeV/c^. A 
five-constraint (5C) fit with an 77 mass constraint is 
used to improve the mass resolution from 8 MeV/c^(4C) 
to 3 MeV/c^, as shown in Fig. [TJc) where xic 40 
is required. To select rj' mesons, \M^+^-j^ — mri>\ < 
0.01 GeV/c^ is required. If more than one combination 
passes the above selection, the combination with A/tt+tt-,; 
closest to m^/ is selected. After the above selection, 
structures similar to those seen for the rj' — >■ 7p chan- 




M...,(GeV/c^) 



M,.,-,(GeV/cO 



uj{u! — )■ 77r")7r"' 

spectively. A clear 77' signal with a 5 MeV/c^ mass resolu- 
tion is evident in the mass spectrum of all selected 77r+7r~ nel in the tt 



FIG. 1: Invariant-mass distributions for the selected candi- 
date events, (a) and (6) are the 77r+7r^ invariant-mass spec- 
trum and the iT^-K~ri' invariant-mass spectrum for rj' — >■ 7p, 
respectively, (c) and (d) are the -k^-k't] invariant-mass spec- 
trum and the -K^-K^-q' invariant-mass spectrum for rj' — >■ 
■K^Tx^T], respectively. The histograms in (b) and (d) are from 
J /ip — >■ 77r"'"7r~77' phase-space MC events (with arbitrary nor- 
malization) for rj' — s- 7p and 77' — s- Tx^-K~r], respectively. 



TT 77' 



invariant mass spectrum can be seen 
in Fig. [IJd), namely peaks near 1.8, 2.1 and 2.4 GeV/c^ 
as well as the 77c. 

Potential background processes are studied with an in- 
clusive sample of 2 x 10® J/ip events generated according 
to the Lund-Charm model Il6i and the Particle Data 
Group (PDG) decay tables [l^- There are no peaking 
backgrounds at the positions of the three resonances. To 
ensure further that the three peaks are not due to back- 
ground, we have studied potential exclusive background 
processes using data. The main background channel is 
from J/t/) ^ tt'' 77+77^77'. Nou-t;' processes are studied 
with 77' mass-sideband events. Neither of these produce 
peaking structures. 

The 7r+7r~77' invariant mass spectrum for the com- 
bined two ri' decay modes is presented in Fig. [21 Here 
a small peak at the position of the /i(1510) signal is also 
present. Fits to the mass spectra have been made us- 
ing four efficiency-corrected Breit-Wigner functions con- 
volved with a Gaussian mass resolution plus a non- 
resonant n^TT^r]' contribution and background represen- 
tations, where the efficiency for the combined channels 
is obtained from the branching-ratio-weighted average 
of the efficiencies for the two 77' modes. The contribu- 
tion from non-resonant "fn^TT^r]' production is described 



by reconstructed MC-generated J/tp — ?> 77r^ 



TT 77 



Phase 



Space (PS) decays, and it is treated as an incoherent pro- 
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cess. The background contribution can be divided into 
two different components, the contribution from non-77' 
events estimated from rj' mass sideband, and the con- 
tribution from J/ip — ;> 7r°7r+7r~77'. For the second back- 
ground, we obtain the background tt^tt^t]' mass spec- 
trum from data by selecting J/ip — > 7r'^7r"'"7r~?7' events and 
reweighting their mass spectrum with a weight equal to 
the MC efficiency ratio of the -fTr^TT^rj' and TT^Tr^Tr^rj' 
selections for J/ip — > 7r°7r+7r^77'. The masses, widths 
and number of event of the /i(1510), the X{1835) and 
the resonances near 2.1 and 2.4 GeV/c^, the X(2120) 
and X(2370), are listed in Table HI The statistical signif- 
icance is determined from the change in — 21nL in the 
fits to mass spectra with and without signal assump- 
tion while considering the change of degree of freedom 
of the fits. With the systematic uncertainties in the 
fit taken into account, the statistical significance of the 
X(1835) is larger than 20a, while those for the /i(1510), 
the X(2120) and the X{2370) are larger than 5.7cr, 7.2cr 
and 6.4(7, respectively. The mass and width from the 
fit of the /i(1510) are consistent with PDG values p7| . 
With MC-determined selection efficiencies of 16.0% and 
11.3% for the rj' jp and rj' Tr'^7r~ rj decay modes 
respectively, the branching fraction for the X{1835) is 
measured to be B(J/i/' ~^ jX{1835)) ■ B{X{1835) 
TT+TT-r]') = (2.87 ± 0.09) x IQ-^. The consistency be- 
tween the two 77' decay modes is checked by fitting their 
7r"'"7r~7y' mass distribution separately with the procedure 
described above. 



TABLE I: Fit results with four resonances for the combined 
two rj' decay modes 




1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 

M,.,-,(GeV/c2) 




1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 

M,.,-,(GeV/c2) 

FIG. 2: (a) The ■k'^-k't]' invariant-mass distribution for the 
selected events from the two 77' decay modes, (b) mass spec- 
trum fitting with four resonances, here, the dash-dot line is 
contributions of non-r?' events and the Tv'^iv'^n~ri' background 
for two 1]' decay modes and the dash line is contributions of 
the total background and non-resonant 7r"'"7r~77' process. 



resonance 


M( MeV/c^) 


r( MeV/c^) 


event 


/i(1510) 


1522.7 ±5.0 


48 ± 11 


230 ± 37 


X(1835) 


1836.5 ± 3.0 


190.1 ±9.0 


4265 ± 131 


X(2120) 


2122.4 ± 6.7 


83 ±16 


647 ± 103 


X(2370) 


2376.3 ± 8.7 


83 ±17 


565 ± 105 



For radiative J/ip decays to a pseudoscalar meson, the 
polar angle of the photon in the J/ip center of mass sys- 
tem, 9^, should be distributed according to 1 + cos^ 9^. 
We divide the jcos^..^! distribution into 10 bins in the 
region of [0, 1.0]. With the same procedure as described 
above, the number of the X(1835) events in each bin 
can be obtained by fitting the mass spectrum in this 
bin, and then the background-subtracted, acceptance- 
corrected I cos 0^1 distribution for the X(1835) is ob- 
tained as shown in Fig. [31 where the errors are statistical 
only. It agrees with 1 + cos^ 9.y, which is expected for a 
pseudoscalar, with x^/d.o.f = 11.8/9. 

The systematic uncertainties on the mass and width 
are mainly from the uncertainty of background represen- 
tation, the mass range included in the fit, different shapes 
for background contributions and the non-resonant pro- 
cess and contributions of possible additional resonances 
in the 1.6 GeV/c^ and 2.6 GeV/c^ mass regions. From 




FIG. 3: The background-subtracted, acceptance-corrected 
|cos6'-,| distribution of the X(1835) for two rj' decay modes 
for J/ip 77r^7r~77'. 



the study of J/ip — ?> pp7r+7r~, the PID efficiency differ- 
ence between data and MC is determined. Using this 
difference and reweighting each MC event with a weight 
equal to the efficiency ratio between data and MC, we 
re-fit the mass spectra and take the changes as system- 
atic uncertainties associated with data and MC incon- 
sistencies for PID efficiencies. The total systematic er- 



rors on the mass and width are 



+5.6 
-2.1 



and 



h38 
-36 



MeV/c 
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for the X(1835), and tfi MeV/c^ for the X(2120), 
14:3 and tf MeV/c^ for the X(2370) respectively. For 
the systematic error of the branching fraction measure- 
ment, we additionaUy inchide the uncertainties of the 
MC generator, charged track detection efficiency, photon 
detection efficiency, kinematic fit, the 77' decay branch- 
ing fractions to 7r+7r~7y and 7p [iTf, the requirement on 
the 77 invariant-mass distribution, signals selection of p, 
77 and rj' and the total number of J/V' events [l3| . The 
main contribution also comes from the uncertainty in the 
background estimation, and the total relative systematic 
error on the product branching fraction for the X(1835) 
is +1^^" 

-18%- 

In summary, the decay channel J/if) — > tt^tt^t]' is an- 
alyzed using two T]' decay modes, ry' — )■ 7/9 and 77' — > 
7r+7r-?7. The X(1835), which was first observed at BESII, 
has been confirmed with a statistical significance larger 
than 20cr. Meanwhile, two resonances, the X(2120) and 
the X(2370) are observed with statistical significances 
larger than 7 .2a and 6.4cr respectively. The masses and 
widths are measured to be: 

• X(1835) 

M = 1836.5 ± 3. 0(stat)t^-j(syst) MeV/c^ 
r = 190±9(stat)j:^^(syst) MeV/c^ 

• X(2120) 

M = 2122.4 ±6.7(stat)l|;^(syst) MeV/c^ 
r = 83 ± 16(stat)j:?i(syst) MeV/c^ 

• X(2370) 

M = 2376.3 ± 8. 7(stat)t'^-j(syst) MeV/c^ 
r = 83 ± 17(stat)l^'^(syst) MeV/c^ 

For the X(1835), the product branching fraction is 
B{J/i} 7X(1835)) • S(X(1835) ^ vr+Tr-Ty') = (2.87 ± 
0.09(stat)^Q'52(syst)) x 10~^, and the angular distribu- 
tion of the radiative photon is consistent with a pseu- 
doscalar assignment. The mass of the X(1835) is con- 
sistent with the BESII result, but the width is signifi- 



cantly larger. If we fit the mass spectrum with one reso- 
nance as BESII, the mass and width of the X(1835) are 
1841.2 ± 2.9 MeV/c^ and 109 ± 11 MeV/c^, where the 
errors are statistical only. 

In the mass spectrum fitting in Fig. [ll^b), possible 
interferences among different resonances and the non- 
resonant process are not taken into account which might 
be a source of the large value for the fit {x^/d.o.f = 
144/62). The dips around 2.2 GeV/c^ and 2.5 GeV/c^ 
may not be fitted well because of the neglect of such 
interferences. In the absence of knowledge of the spin- 
parities of the resonances and their decay intermediate 
states, reliable fits that include interference cannot be 
done. To determine the spin and parity of the X(1835), 
X(2120) and X(2370), and to measure their masses and 
widths more precisely, a partial wave analysis must be 
performed, which will be possible with the much higher 
statistics J/ip data samples planned for future runs of 
the BESIII experiment. 
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